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Abstract 



The multidimensional braneworld gravity model by Dvali, Gabadadze 
and Porrati was primarily put forth to explain the observed accelera- 
tion of the expansion of the Universe without resorting to dark energy. 
If^ \ One of the most intriguing features of such a model is that it also pre- 

■ diets small effects on the orbital motion of test particles which could 
^ I be tested in such a way that local measurements at Solar System scales 
O^' would allow to get information on the global properties of the Universe. 
5^ I Lue and Starkman derived a secular extra-perihelion uj precession of 

5 X 10^'' arcseconds per century, while lorio showed that the mean 
longitude A is affected by a secular precession of about 10"'^ arcsec- 
, . ■ onds per century. Such effects depend only on the eccentricities e of 

r> I the orbits via second-order terms: they are, instead, independent of 

■ their semimajor axes a. Up to now, the observational efforts focused 

on the dynamics of the inner planets of the Solar System whose orbits 
are the best known via radar ranging. Since the competing Newtonian 
and Einsteinian effects like the precessions due to the solar quadrupole 
mass moment J2, the gravitoelectric and gravitomagnetic part of the 
equations of motion reduce with increasing distances, it would be pos- 
sible to argue that an analysis of the orbital dynamics of the outer 
planets of the Solar System, with particular emphasis on Saturn be- 
cause of the ongoing Cassini mission with its precision ranging instru- 
mentation, could be helpful in evidencing the predicted new features 
of motion. In this paper we investigate this possibility by comparing 
both analytical and numerical calculation with the latest results in the 
planetary ephemeris field. Unfortunately, the current level of accuracy 
rules out this appealing possibility. 
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1 Introduction 

Recently, Dvali, Gabadadze and Porrati (DGP) put forth a model of grav- 
ity ^ which allows to explain the observed accelerated expansion of our 
Universe without resorting to the concept of dark energy. In such a picture 
our Universe is a (3+1) space-time brane embeddeed in a 5-dimensional 
Minkowskian bulk. While the dynamics of the Standard Model particles 
and fields is confined to our brane, gravity can fully explore the entire bulk 
getting strongly modified at large distances of the order of ro ~ 5 Gpc. 
An intermediate regime is set by the Vainshtein scale = {rgr^)^^^, where 
rg = 2GM /(? is the Schwarzschild radius of a central object acting as source 
of gravitational field. For a Sun-like star amounts to about 100 parsec. 
An updated overview of the phenomenology of DGP gravity can be found 
in 12]. 

1.1 Local orbital effects 

One of the most interesting features of such a picture is that in the process 
of recovering the 4-dimensional Newton-Einstein gravity for r « « tq, 
DGP predicts small deviations from it which yield to effects observable at 
local scales l3l. They come from an extra radial acceleration of the form 



The minus sign is related to a cosmological phase in which, in absence of 
cosmological constant on the brane, the Universe decelerates at late times, 
the Hubble parameter H tending to zero as the matter dissolves on the brane: 
it is called Friedmann-Lemaitre-Robertson- Walker (FLRW) branch. The 
plus sign is related to a cosmological phase in which the Universe undergoes 
a deSitter expansion with the Hubble parameter H = c/ro even in absence of 
matter. This is the self-accelerated branch, where the accelerated expansion 
of the Universe is realized without introducing a cosmological constant on 
the brane. Thus, there is a very important connection between local and 
cosmological features of gravity in the DGP model. About the local effects, 
Lue and Starkman |S] derived an extra-secular precession of the pericentre ui 




(1) 
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of a nearly circular orbit of a test particle of 5 x 10^^ arcseconds per century 
(" cy~^), while lorio [H] showed that also the mean anomaly Ai is affected 
by DGP gravity at a larger extent; the longitude of the ascending node is 
left unchanged. As a result, the mean longitude X = uj + Q + M, which is 
a widely used orbital parameter for nearly equatorial and circular orbits as 
those of the Solar System planets, undergoes a secular precession of the order 
of 10~^ " cy~^. It is independent of the semi-major axis o of the planetary 
orbits and depends only on their eccentricities e via second-order terms. 
Recent improvements in the accuracy of the data reduction process for the 
inner planets of the Solar System [71 |HI, which can be tracked via radar- 
ranging, have made the possibility of testing DGP very thrilling [Qj l6l ITn]. 
In particular, lorio |10J showed that the recently observed secular increase 
of the Astronomical Unit jlH I12j can be explained by the self-accelerated 
branch of DGP and that the predicted values of the Lue-Starkman perihelion 
precessions for the self-accelerated branch are compatible with the recently 
determined extra-perihelion advances 0, especially for Mars, although the 
errors are still large. 

1.2 Problems of DGP gravity 

There are still some issues about the theoretical consistency of the DGP 
model, so that direct, local observational tests would be very useful and im- 
portant to tackle observationally this matter. The recovery of the Newton- 
Einstein four dimensional gravity is rather not trivial. The aforementioned 
orbital effects come from the correction to the Newtonian potential of a 
Schwarzschild source found in [11121. Such a potential is obtained within a 
certain approximation which is valid below the Vainshtein scale. However, 
it is not yet clear, at present, whether this potential can match contin- 
uously onto a four-dimensional Newtonian potential above the Vainshtein 
scale, and then, also match onto the five-dimensional potential above the 
crossover scale tq. An alternative solution that smoothly interpolates be- 
tween the different regions was discussed in I14j . The correction to the 
Newtonian potential arising from that solution below the Vainshtein scale 
is somewhat different from what used here. In particular, it is reduced by a 
multiplicative factor smaller than unity. As a consequence, the predictions 
are also different. Moreover, according to the authors of JSI, ghosts affect 
the self-accelerated branch, plaguing the consistency of the model at cosmo- 
logical scales. In regard to the domain below the Vainshtein scale, according 
to a ghost appears in the Pauli-Fierz model of massive gravity, while it 
is absent in DGP on the conventional branch. A similar result has also been 
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obtained in [T7| . 

1.3 Aim of the paper 

As already pointed out, up to now the only available preliminary tests of the 
DGP gravity refer to the inner planets of the Solar System. In this paper 
we focus on the outer planets in order to enlarge the possible set of inde- 
pendent checks. For such celestial bodies many competing Newtonian and 
Einsteinian effects are smaller than the DGP features of motion allowing, in 
principle, cleaner tests of the Lue-Starkman and lorio precessions. 

2 Can the outer planets of the Solar System be 
useful? 

2.1 The DGP precessions and the competing Newtonian and 
Einsteinian effects 

We list the DGP precessions and the competing Newtonian and Einsteinian 
precessions in Tabled taken from p|, for the mean longitudes. In order 
to extract the DGP signal from the analysis of the determined secular pre- 
cessions of the orbital elements it is mandatory that the competing effects 
are included in the dynamical force models of the planetary data reduction 
softwares to a sufficient level of accuracy. It occurs for all such features of 
motion apart from the Lense-Thirring effect, which is not modelled at all 
in the currently used ephemeris softwares, and the solar quadrupole mass 
moment J2 which is, in fact, modelled but whose uncertainty is of the order 
of 10%. Such dynamical features induce residual mismodelled precessions 
which could seriously affect the recovery of a smaller effect. Otherwise, it 
could be possible to suitably combine the orbital elements in order to a priori 
cancel out some of the unwanted perturbations. 

A careful inspection of Table ^ shows that for the outer planets of the 
Solar System-and for Mars-the nominal values of the LT and J2 precessions 
are themselves well smaller than the DGP rate. The nominal GE preces- 
sions are larger but by a sufficiently small amount to create no problems 
because they are routinely included in the dynamical force models in terms 
of the PPN parameters f3 and 7 |^ which are known with an accuracy of 
10~^ ^ or better ^23; • Moreover, deviations from their relativistic values 
are expected at the level of 10~^ — 10~^ "^^^ Newtonian N-body pre- 
cessions are known very well and the impact of the asteroids, which limits 
the obtainable accuracy for the inner planets with particular emphasis on 
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Table 1: Nominal values, in " cy~ , of the secular precessions in- 
duced on the planetary mean longitudes A by the DGP gravity and 
by some of the competing Newtonian and Einsteinian gravitational per- 
turbations. For a given planet, the precession labelled with Numeri- 
cal includes all the numerically integrated perturbing effects of the dy- 
namical force models used at JPL for the DE200 ephemeris. E.g., it 
also comprises the classical N- body interactions, including the Keple- 
rian mean motion n. For the numerically integrated planetary preces- 
sions see http://ssd.jpl.nasa.gOv/elem_planets.html#rates. The effect la- 
belled with GE is due to the post-Newtonian general relativistic gravito- 
electric Schwarzschild component of the solar gravitational field TS;, that 
labelled with J2 is due to the Newtonian effect of the Sun's quadrupole mass 
moment J2 |18j and that labelled with LT is due to the post-Newtonian gen- 
eral relativistic gravitomagnetic Lense-Thirring [20] component of the solar 
gravitational field (not included in the force models adopted by JPL). For 
J2 the value 1.9 x 10"'' has been adopted p. For the Sun's proper angular 
momentum S, which is the source of the gravitomagnetic field, the value 



1.9 X 10^1 


kg m^ s 


1 


1211 has been adopted. 










Planet 


DGP 




Numerical 


GE 


J2 


LT 




Mercury 


1.2 X 10" 


-3 


5.381016282 x 10^ 


-8.48 x 10^ 


4.7 X 10-2 


-2 X 10" 


-3 


Venus 


1.2 X 10" 


-3 


2.106641360 x 10^ 


-1.72 x 10^ 


5 X 10-3 


-3 X 10" 


-4 


Earth 


1.2 X 10" 


-3 


1.295977406 x 10^ 


-7.6 


1.6 X 10-3 


-1 X 10" 


-4 


Mars 


1.2 X 10" 


-3 


6.89051037 x lO'^ 


-2.6 


3 X 10-^ 


-3 X 10" 


-5 


Jupiter 


1.2 X 10" 


-3 


1.09250783 x lO"^ 


-1 x 10-^ 


5 X 10"*^ 


-7 X 10" 


-7 


Saturn 


1.2 X 10" 


-3 


4.4010529 X 10^ 


-2 x 10-2 


6 X 10-^ 


-1 X 10" 


-7 


Uranus 


1.2 X 10" 


-3 


1.5425477 x 10^ 


-4 x 10-3 


5 X 10-*^ 


-1 X 10" 


-8 


Neptune 


1.2 X 10" 


-3 


7.864492 x 10^ 


-1 x 10-3 


1 X 10-s 


-5 X 10" 


-9 


Pluto 


1.2 X 10" 


-3 


5.227479 x 10^ 


-7 X 10-^ 


4 X 10-^ 


-2 X 10" 


-9 
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Table 2: Formal standard statistical errors in the non-singular planetary 
orbital elements, from Table 4 of [7j- The latest EPM2004 ephemerides 
have been used. The angles i and w are the inclination and the longitude 
of perihelion, respectively. The realistic errors may be up to one order of 
magnitude larger. The units for the angular parameters are milliarcseconds 



(mas). 


Planet 


a (m) 


sin i cos (mas) sin i sin (mas) e cos w 


(mas) esinro (mas) 


A (mas) 


Mercury 


0.105 


1.654 


1.525 


0.123 


0.099 


0.375 


Venus 


0.329 


0.567 


0.567 


0.041 


0.043 


0.187 


Earth 


0.146 






0.001 


0.001 




Mars 


0.657 


0.003 


0.004 


0.001 


0.001 


0.003 


Jupiter 


639 


2.410 


2.207 


1.280 


1.170 


1.109 


Saturn 


4222 


3.237 


4.085 


3.858 


2.975 


3.474 


Uranus 


38484 


4.072 


6.143 


4.896 


3.361 


8.818 


Neptune 


478532 


4.214 


8.600 


14.066 


18.687 


35.163 


Pluto 


3463309 


6.899 


14.940 


82.888 


36.700 


79.089 


Mars, is. 


of course. 


negligible. So, 


at first sig 


;ht, it would seem 


appealing 





to consider such planets as ideal candidates to test DGP gravity because of 
the absence of competing aliasing effects which could mask it. This feeling 
is enforced by the ongoing Cassini tour of the Saturn system which started 
at the end of 2004 and should last for almost five years. Lue 2 hoped that 
the precise radio-link to Cassini might, indeed, provide ideal test for the 
anomalous periheion precession. 

2.2 Confrontation with data 

Unfortunately, the situation is in fact rather unfavorable for the outer plan- 
ets. Indeed, it must be pointed out that the investigated effect is a secular 
one, i.e. it is integrated over one orbital period of the planet. The currently 
available modern observations for the outer planets span almost one century 
and cover a limited number of full orbital revolutions of Jupiter (P = 11.86 
yr), Saturn {P = 29.46 yr) and Uranus (P = 84 yr). Up to now, the orbit 
of Jupiter is the best known among the outer planets because a number 
of precise radar observations by spacecraft (Pioneer 10 and 11, Voyager 1 
and 2, Ulysses and Galileo) approaching the planet or orbiting it have been 
performed. The other planets rely entirely upon optical observations. The 
latest results by Pitjeva jTj are reported in Table HI It can be noted that, 
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even in the case of Jupiter, the currently available sensitivity does not al- 
low to detect the DGP precession. In regard to the mean longitudes of 
Jupiter and Saturn, their (formal) uncertainties are slightly larger than the 
lorio precession. Another source of limiting systematic bias is represented 
by the uncertainty 6n = (3/2) Y^GM/a^(5a in the Keplerian mean motion 
n = \J GMjcfi induced by the errors ba in the semimajor axis. According to 
Table 121 for Jupiter and Saturn we have bn ~ 10~^ " cy~^ which is one order 
of magnitude larger than the investigated effect. The situation is worse for 
the perihelia: indeed, the eccentricities of Jupiter and Saturn only amount 
to 0.04 and 0.05, respectively, so that their apsidal lines are very difficult 
to be determined. As a consequence, from Table |21 it can be obtained an 
uncertainty up to 100 times larger than the Lue-Starkman perihelion preces- 
sion. It is important to stress that such evaluations are based on the formal, 
statistical errors in the planetary orbital elements: realistic errors may be 
up to ten times larger. 

2.2.1 Numerical simulations 

The Keplerian orbital elements are not directly measured quantities: they 
are related to data in an indirect way. For the outer planets the true observ- 
ables are the right ascension a and the declination b. Figure 2 of ^ shows 
the residuals, in of a cos 5 and b for Jupiter, Saturn, Uranus, Neptune 
and Pluto; the scale is ±5 ". In order to yield a direct and unambiguous 
confrontation with the observations we decided to numerically produce a 
set of residuals by integrating the planetary equations of motion with and 
without the DGP acceleration of eq.(p3) with the plus sign. 

To integrate the equations of motion the Mercury package |25j has been 
chosen among the software packages freely available on the Web because it 
is of simplest use and easily allows to introduce user-defined perturbative 
forces, as the DGP acceleration of eq.(p3). This package has been also used 
in EH]. 

The Mercury package can use the following integration methods 

• second order mixed-variable symplectic 

• Bulirsh-Stoer (general) 

• Bulirsh-Stoer (conservative systems) 

• Everhardt Radau 15th order 

• hybrid (symplectic / Bulirsh-Stoer) 
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The perturbative force is introduced by giving the x, y, and z components 
in an heUocentric frame of the corresponding acceleration. 

Thanks to the experience gained during the integrations for yi'6\ , the Ev- 
erhardt Radau method has been chosen as reference, due to its precision and 
speed. The initial epoch is 1913.0 (JD 2419770.5), and the initial planetary 
positions have been obtained from JPL ephemerides. 

The integration has been performed twice, with and without the DGP 
acceleration of eq.(^, and a cos 5 and 5 have been calculated in both cases. 
The difference between a cos 6 and 6 with and without the perturbation 
constitute the final results. They are represented in Figure ITl-Figure 1101 

As a control, the same procedure has been performed adopting the 
Bulirsh-Stoer (general) method of integration; the results are not reported 
here, because the resulting figures are exactly equal to those here presented 
(the differences are at level of the sixth significative digit). 

As can be noted, the pattern which would be induced by a DGP pertur- 
bation on the planetary motions cannot be discerned with the present-day 
orbital accuracy. 

2.3 The role of Cassini and of other future interplanetary 
missions 

In regard to Cassini, it turns out that recent attempts performed at JPL to 
use its first data to improve the orbit of Saturn failed to reach the required 
level of accuracy (Jacobson, R.A., private communication, 2005). Moreover, 
the duration of the Cassini mission will only cover one sixth of one entire 
orbital period of Saturn. In regard to^ GAIA, it is not tailored for observa- 
tions of planets: indeed, they are too bright and cannot be measured exactly. 
Moreover, as in the case of Cassini, long stream of observations should be 
required, not mere sparse points. Thus, it is doubtful that any accurate 
observations of planets will be obtained from GAIA (Pitjeva, E.V., private 
communication, 2005). 

As a consequence, the preferred test-bed for DGP will likely remain 
the inner Solar System even in the near future when the data from the^ 
BepiColombo, Messenger, Venus Express spacecraft will be available and 
especially if the currently investigated projects on the interplanetary laser 
ranging will be implemented. Indeed, the use of optical frequencies 

^See on the WEB http://sci.esa.mt/science-e/www/area/mdex. cfm?fareaid=26l 
^See on the WEB http://sci.esa. int/science-e/www/area/index.cfm?fareaM=30] 
for BepiColombo, http://messenger.jhuapl.edu/ for Messenger and 
http://sci.esa. int/science-e/www/area/index.cfm?fareaid=64 for Venus Express. 
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Jupiter 
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Figure 1: Jupiter: DGP shift, in arcseconds, for a cos (5 over T = 90 years. 
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Jupiter 
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Figure 2: Jupiter: DGP shift, in arcseconds, for S over T = 90 years. 
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Saturn 
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Figure 3: Saturn: DGP shift, in arcseconds, for a cos (5 over T = 90 years. 
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Figure 4: Saturn: DGP shift, in arcseconds, for S over T = 90 years. 
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0.02 



Uranus 




-0.03 I 1 1 1 1 1 

1913 1931.2 1949.4 1967.6 1985.8 2004 

Time (yr) 



Figure 5: Uranus: DGP shift, in arcseconds, for a cos 5 over T = 90 years. 
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Figure 6: Uranus: DGP shift, in axcseconds, for S over T = 90 years. 
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Figure 7: Neptune: DGP shift, in arcseconds, for a cos S over T = 90 years. 



15 




Figure 8: Neptune: DGP shift, in axcseconds, for 6 over T = 90 years. 
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Figure 9: Pluto: DGP shift, in axcseconds, for a cos S over T = 90 years. 
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Figure 10: Pluto: DGP shift, in axcseconds, for 6 over T = 90 years. 
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would definitely allow to overcome the limitations posed by the solar corona 
to the radar ranging technique. 

3 Conclusions 

In this paper we analyzed the possibility of using the outer planets of the 
Solar System, especially Jupiter and Saturn, to test the Dvali-Gabadadze- 
Porrati multidimensional braneworld model of gravity via the Lue-Starkman 
secular precession of perihelion u; and the lorio secular precession of the mean 
longitude A. Indeed, for the distant planets such effects, which are indepen- 
dent of the geometrical features of the orbits apart from second-order terms 
in the eccentricities e, would be larger than the other competing perturba- 
tions induced by the mismodelled or unmodelled Newtonian and Einsteinian 
features of motion. Moreover, the current presence of the Cassini space- 
craft in the Saturnian system with its precise radiotcchnical apparatus and 
the forthcoming astromctric mission GAIA might induce some expectations 
about such an intriguing possibility. Unfortunately, the real situation is less 
favorable than it was hoped. Indeed, the investigated new features of mo- 
tion are currently by one- two orders of magnitTidc below the threshold set 
by the (formal) accuracy of the most recent determinations of the orbital 
elements of Jupiter and Saturn. Moreover, it is doubtful that Cassini will 
substantially contribute to improving our knowledge of the orbit of Saturn 
to a sufficient level. Indeed, it could yield only sparse points because its 
lifetime will span just one sixth of an entire orbital period of Saturn. 
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